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Abstract 



C^ . Non-equihbrium phenomena occur not only in physical world, but also in finance. In 

2 I this work, stochastic relaxational dynamics (together with path integrals) is applied 

F^ I to option pricing theory. Equilibrium in financial markets is defined as the absence 

^ • of arbitrage, i.e. profits "for nothing". A recently proposed model (by Ilinski et al.) 

considers fiuctuations around this equilibrium state by introducing a relaxational 
dynamics with random noise for intermediate deviations called "virtual" arbitrage 
^ I returns. In this work, the model is incorporated within a martingale pricing method 

^ I for derivatives on securities (e.g. stocks) in incomplete markets using a mapping 

Q>^ ■ to option pricing theory with stochastic interest rates. The arbitrage return is 

considered as a component of a fictitious short-term interest rate in a virtual world. 
Q . The infiuence of intermediate arbitrage returns on the price of derivatives in the 

G^ i real world can be recovered by performing an average over the (non-observable) 

arbitrage return at the time of pricing. Using a famous result by Merton and with 
c^ ■ some help from the path integral method, exact pricing formulas for European call 

and put options under the infiuence of virtual arbitrage returns (or intermediate 
,-rt ■ deviations from economic equilibrium) are derived where only the final integration 

^ ! over initial arbitrage returns needs to be performed numerically. This result, which 

has not been given previously and is at variance with results stated by Ilinski et al., 
is complemented by a discussion of the hedging strategy associated to a derivative, 
which replicates the final payoff but turns out to be not self-financing in the real 
world, but self-financing when summed over the derivative's remaining life time. 
5^ i Numerical examples are given which underline the fact that an additional positive 

risk premium (with respect to the Black-Scholes values) is found reflecting extra 
hedging costs due to intermediate deviations from economic equilibrium. 



1 Introduction 

The pricing and hedging of derivatives is a major task for financial institutions pS 



and has become an increasingly popular topic in statistical physics 1131 , 0. Deriva- 



tives are sometimes also called contingent claims, as the buyer of the derivatives is 
entitled to receive a certain payoff up to (or at) some future time T, the time of 
expiry, dependant on the price 5* of a so-called "underlying" security (say a stock) 
within a certain time interval between today and T or at time T. The simplest case 
is a so-called European call (put) option which gives the buyer the right to buy (sell) 
an underlying security (e.g. a stock) at a certain time T in the future for a fixed 
price K (the strike price). These options are also called "plain vanilla options" for 
their simplicity (as common as the vanilla icecream flavor). 

The classical result of Black and Scholes [Q on option pricing which revolution- 
ized the world of finance and still forms the foundation for most of modern research, 
is based on the existence of an equilibrium, generally called "absence of arbitrage" , 
i.e. the impossibility of a profit "for nothing" . 

The use of the no-arbitrage assumption for pricing purposes is nicely elucidated in 



a standard text book like [^ where simple pricing equations for forward contracts 
are derived from optimization arguments. If e.g. the forward price F at time t 
for buying or selling (assuming no bid/offer spreads and transactions costs) a non- 
dividend paying security S" at a later time T were less than S'exp(r(T — t)), then 
a riskless profit could be obtained in the following way: at time t, one enters into 
a forward contract to buy the security for F at time T, and one short sells the 
security (i.e. one borrows the security from somebody else and sells it, assuming no 
fees for simplicity) and puts the proceeds on a deposit at the riskless interest rate 
r (assuming no credit risk); at time T, one receives the security from the forward 
contract thereby closing out the short position in the security (i.e. handing it over to 
the lender), while receiving the nominal amount plus interest from the deposit less 
the forward price paid: S'exp(r(T — t)) — F > 0. Likewise, an arbitrage is possible 
when F > Sexp[r{T — t)). As the information on either situation spreads in the 
market, the inequalities disappear, and the relation F = S exp{r{T — t)) results. 
Now the no-arbitrage assumption anticipates the equality to hold right from the 
beginning, thus implying that arbitrage opportunities disappear infinitely fast. Now 
many trading activities are motivated exactly by the fact, that this is not the case, 
but that arbitrage returns exist in the market for a short time Tarbitrage > 0. After 
this time, the information on arbitrage opportunities has reached enough market 
participants to make them disappear. 

The absence of arbitrage assumption paves the way for one of the fundamental 



pillars of mathematical finance which is the theorem by Harrison and Pliska ^\. 
In fact whenever markets are complete (i.e. when any derivative can be hedged 
by a self-financing strategy, which is a more restrictive statement than absence of 
arbitrage), then there is a unique equivalent martingale measure for the underlying 
security and vice versa (see [^ for an introductory discussion on martingale theory). 



A stochastic process Xt is a martingale with respect to the measure Q if and only 

iiEQ[\Xt\] <oo 

X, = Eq[X,|J-,], s>t (1) 

where Tt is the filtration at time t, i.e. the information accumulated until time t. 



This rather technical statement is the basis for risk-neutral valuation: Derivatives 
can be priced in a world where all yields are equal to the risk-free interest rate 
(minus dividend yields etc.). 

Apart from the dynamic deviations from the no-arbitrage situation discussed 
above, there is a large literature on serious drawbacks of the Black-Scholes model 
itself which is classically used to implement no-arbitrage pricing of options, i.e. that 
price returns evolve according to Brownian motion with constant drift and volatility. 
Empirical studies of return distributions in fact show volatility clustering and fat 
tails IQ, D, 1^, Q, and so real price changes appear to be more efficiently modelled 
by truncated Levy processes (TLP) [|, ^j. However, rational option pricing using 
the martingale approach appears to be still working (see ^ and ref.s therein), so 
the main theme of the Black-Scholes method, i.e. the possibility to set up a self- 
financing hedging strategy (a notion to be explained below) seems to hold. Moreover, 
studies on autocorrelations of price changes or on the distributions of price changes 
themselves (taken for different time scales) demonstrate a crossover to gaussian 
dynamics after a certain time scale which might vary from several minutes to days 
(depending on market liquidity) [0, |14 . 

In the present work, this time scale is proposed to be proportional to Tarbitrage- In 
this sense, fat tails of distributions of price changes are a signature of intermediate 
arbitrage opportunities on short time scales (compared to Tarbitrage)- The implica- 
tions for pricing options that are not very short-lived seem to be that intermediate 
arbitrage opportunities may be modelled as deviations from Brownian motion, and 
thus as perturbative effects. One way to treat these deviations are stochastic volatil- 
ity (SV) models (see for a review in the context of option pricing [Q). The present 
work complements these models and gives an alternative approach which modifies 
the drift of the asset price process rather than its volatility. 

If one tries to get rid of the drawbacks of the Black-Scholes model by dynamic 
parameters which are not directly tradable such as in SV models or in the approach 
discussed below, one encounters a new problem: as opposed to complete markets 
defined above, a self-financing strategy using traded instruments ceases to exist. In 
general, any hedging strategy can only reduce the risk inherent in the final payoff to 
an intrinsic component 0, |13|. Technically, this leads to more than one equivalent 



martingale measure |]T^. ambiguity for derivatives pricing is handled by introducing 
additional constraints on the hedging strategy, e.g. minimizing the expected squared 
cost for the remaining life time of the option while exactly replicating the final payoff 
(local risk-minimization) or minimizing the expected squared net loss at the time 
of maturity of the option (mean-variance hedging) 0, H, |16| (see also |jl3|, |15|] for a 
physicist's approach). For the model presented in this work, a very specific method 
is proposed in order to select an equivalent martingale measure which satisfies both 
constraints simultaneously. 

The issue of option pricing in incomplete markets has become a matter of prac- 
tical interest recently, in particular with the increasing importance of credit deriva- 
tives. As opposed to conventional derivatives which cover market risks, there is not 
a large underlying market of actively traded credit risk instruments in every credit 
risk category. As opposed to a stock, e.g. a loan is usually not traded. 

Considering fully complete and incomplete markets as extreme cases of real mar- 
kets, one is naturally forced to ask for crossover effects or transitions between the 



two regimes. A possible answer to this question might be given in terms of a dy- 
namic model which considers market incompleteness in terms of fluctuations around 
an economic equilibrium characterizing a complete market. 

An important step in this direction has been given by Ilinski and Stepanenko |T0[ 



and Ilinski ||TT| who assume the existence of intermediate, "virtual" arbitrage returns 
Xt, which appear and disappear over a certain time scale which may be identified 
with Tarbitrage mentioned above. In fact, Ilinski et al. intend to treat arbitrage 
effects as a perturbation to the usual Black-Scholes risk-free rate r, which gives the 
yield on all investments in a risk-neutral world. More specifically, the Black-Scholes 
risk-free rate is split into a constant part r" and an arbitrage return xt according to 
ft = r^ + Xt- This model is different from stochastic volatility models: it modifies the 
(risk-adjusted) drift of the asset price process rather than its volatility. In principle 
there appears no reason to favor one approach over the other. In fact, preliminary 
results from simulations of the stochastic drift model versus the stochastic volatility 
model by Stein & Stein |]T^ seem to point in this direction. Detailed results on this 



comparison will be published elsewhere [^. The stochastic drift approach discussed 



here has the advantage that it can be mapped to models with stochastic interest 
rates (see below). 



Ihnski et al. [|T3l[|III present two approaches to calculate option prices. In [|T0 
a perturbative method is given based on the classical Black-Scholes equation where 
the constant risk-free rate r is replaced by r° + xt where r° is constant and Xt is 
random. This equation is then iterated to second order in xt and averaged over 
with respect to Xt. The results obtained for European call and put options are not 
reproduced by our exact calculation. We come back to this difference at the end of 
section 5 which explains why our results are reasonable for the specific dynamics of 
the arbitrage return used here. The observed difference just mentioned leads us to 



suggest that either the method or some further approximations that made in |10 



are not correct. In a second paper |TT[ Ilinski et al. derive a fully deterministic 
Black-Scholes type PDE that depends both on the current level of the security price 
St and the arbitrage return xt. This equation is not further evaluated. In section 3, 
we show that it contains a misprint. We add as a remark that the issue of measure 
change and the construction and selection of an equivalent martingale measure which 
is fundamental to option pricing (see e.g. [^) is not addressed. 

Therefore, in this article a different route is proposed. First, the arbitrage return 
Xt is considered as a part of a stochastic interest rate dynamics for the risk-free rate 
rt in a virtual world (where arbitrage returns are directly observable). Essentially, 
one is led to a Black-Scholes type equation for a derivative depending on two state 
variables, the security price St and the arbitrage return xt which is derived from the 
risk-free rate in the virtual world according to rt = r^ + Xt. The constant part r^ is 
supposed to be the risk-free rate in the real world which consequently is assumed to 
be constant. The reason for the latter simplification is the later comparison with the 
Black-Scholes pricing formulas. The implementation of Xt as a part of a fictitious 
interest rate process leads to a stochastic drift for the asset process St (with respect 
to the particular martingale measure chosen, for details see below) and thus couples 
the dynamics of Xt to 5*^ As the arbitrage return is an intermediate phenomenon 
on time scales shorter than the time to expiry, we follow Ilinski et al. by enforcing 
the boundary condition at the time of expiry of the option that the arbitrage return 



should disappear. This constraint may be relaxed, however, if one allows for the 
possibility that any hedging strategy might not replicate (i.e. provide for) the final 
payoff of the option. Nonetheless, we stay with this constraint (also in order to 
compare with the results of Ilinski et al. [|TO||]TT|). However, we do not implement 



this condition into the payoff function of the option like Ilinski [|ri|], but into the 



average over arbitrage returns. This procedure allows us to use a famous result by 



Merton on options in a stochastic interest rate environment |26|. It is not meant to 
imply that intermediate arbitrage returns can be thought of as the random part of 
real interest rates. After averaging, we obtain a previously unknown exact result for 
European claims under the influence of virtual arbitrage. (The pricing of American 
claims which may be exercised prior to their maturity is possible in principle using 
a "tree" procedure [^, i.e. a scheme based on discrete probabilities and discrete 
time). These exact pricing formulas differ significantly from the results obtained by 
Ilinski et al. |1^ . 



The outline of this paper is as follows: In the next section, we present the route 
from the Black-Scholes model to a non-equilibrium market model, taking up the idea 
of intermediate ( "virtual" ) arbitrage by Ilinski and Stepanenko. The third section 
will show how the effect of arbitrage returns on option pricing can be considered 
in terms of a stochastic interest rate environment in a virtual world. In section 4, 
European call and put options are valued in the presence of virtual arbitrage returns. 
In section 5, the issue of a replicating hedging strategy both in the virtual and real 
world and the selection of an equivalent martingale measure is addressed. Some 
explicit numerical pricing examples are given and their difference to the classical 
Black-Scholes results are explained in section 6. In the final section, the results are 
briefly discussed. 

2 From the Black-Scholes model to the dynamics 
of arbitrage returns 

Let us briefly review the Black-Scholes analysis in order to motivate the notion of 
arbitrage returns, following |TI]|. The model for a one security market is given by 



dSt = i^Stdt + aStdWl (2) 

where St is the security price, /i the drift, and dW^ a Wiener process. It may be 
motivated from the fact that ln(S'j+i/5'j), where i + 1 and i denote discrete points 
in time, performs a random walk [Q. Now the price of a derivative Vj = V{St,t) 
whose payoff is contingent on the security price St at some future time T can be 
determined by setting up a portfolio fl^ consisting of the derivative Vt and a position 
—A of the security St'. 

nt = Vt- ASt (3) 

Then if A = |^, where S = St, this portfolio is riskless as uncertainties arising from 
the Wiener process are eliminated which can be seen by evaluating dUt using Ito's 
lemma. Therefore, the portfolio is known to grow at the risk-free rate, i.e. 

dUt = rUtdt (4) 



For constant interest rates, equating expressions for dUt gives the Black-Scholes 
partial differential equation (PDE): 

Specifying a certain boundary condition to this equation representing the option 
payoff at the time of maturity completes the usual Black-Scholes pricing problem. 
As a reminder, let us note that the drift fi which was introduced in the market model 
Eq.(ll) is absent from the pricing equation Eq.(|^). 

The idea of arbitrage returns may be motivated by assuming that in the presence 
of arbitrage opportunities, the true return of the portfolio Ut is not equal to the 
constant risk-free interest rate r, but might be less or more than that. Following 



TT[| , an arbitrage return xt is now introduced by substituting for r 

Tt = r^ + Xt (6) 

where Xt is assumed to follow the dynamics of a decay process with random noise: 

-Xxt + Tjt (7) 



dxt 



dt 
where rjt is characterized by: 

(r?t) = 0, {'ntr„)=T.H{t-t') (8) 



As to the nature of rjt further complications are discussed in |]T0| , but they are not 
important for our analysis. Basically, a stochastic component Xt as been added to 
the constant risk-free rate tq. The question now is: How does the process for the 
arbitrage return xt affect the price of a derivative ? 

Substituting Eq.(|^) for risk-free rate r in the standard Black-Scholes PDE, Ilinski 
and Stepanenko simply proceed and derive the following PDE: 



CbsV = xt\V- S— (9) 




where Cbs is the operator from the standard Black-Scholes PDE, CbsV = 0, for 
r = r°. We will clarify below that the replacement r ^ r° + x^ in fact is equivalent 
to introducing an interest rate process rt = r^ + Xt in a virtual world where tradable 
instruments dependant on this interest rate exist. 

The specific origin of intermediate arbitrage returns and market incompleteness 
is assumed to be contained in the parameters A which sets the time scale for de- 
viations from market equilibrium and S which gives a measure for the arbitrage 
returns themselves. Specific values are discussed in section 6. Of course, in reality 
transaction costs which are neglected here for simplicity might effectively destroy 
small arbitrage returns. 

In the next section, a different approach is presented in the framework of standard 
option pricing theory which allows to study the influence of intermediate deviations 
from financial equilibrium (as defined by the no-arbitrage assumption) on derivative 
pricing. 



3 Derivatives in the presence of arbitrage oppor- 
tunities: a mapping to option pricing theory 
with stochastic interest rates 

The way the arbitrage return Xt has been introduced in the last section, in particular 
that the portfolio Ut grows at the rate r^ + Xt, allows for a mapping to option pricing 
theory with stochastic interest rates. We will call rt = r^ + Xt an interest rate in 
a virtual world, but do not insinuate that the arbitrage return is a part of the real 
interest rate. This virtual world will serve as a stage where known results can be 
used, but finally these results need to be projected to the real world. Let us for the 
moment assume, that this virtual world can be set up. The justification for its use 
will be delayed to section 5. 

Let us first review the PDE approach to option pricing with stochastic interest 
rates. The stochastic nature of a (short) interest rate vt is usually taken into account 
by stating a stochastic differential equation (SDE) as follows: 

dn = pin, t)dt + S(n, t)dw^ (10) 

The parameters p and S specify drift and volatihty respectively, and may depend on 
rt and t. The drift specifies the deterministic ( "trend" ) component of the interest 
rate dynamics whereas the volatility describes the stochastic fluctuations. The in- 
crement dWf is a Wiener process. The general PDE for a derivative V = V{S, r, t) 
dependant on S = St and r = r^ can be found in the literature p3|. Assuming for 
simplicity no correlations between the Wiener processes dW^ from Eq.(^ and dW^ 
and suppressing the functional dependance of S and p, the PDE is given by: 

dv a's^a^-v ^dv, , , „ Eav av, , „ „ 

IF + — as? ^ ^as <" " ^■'^> -'""^ll^^ * <" " ^^^' = » <"' 

The parameters Aj, i = 1,2 are known as the market prices of risk for the security 
S and the risk- free rate r. They can be obtained by finding the change of measure 
which makes the respective discounted price process a martingale ^^. For a non- 



dividend paying security governed by Eq.(|D, Ai = (/x — r)/a. Incidentally, if r is 
constant, one recovers the Black-Scholes PDE Eq.(P). As r is not a tradable security, 
a tradable interest rate instrument is needed, e.g. a zero bond with maturity T whose 
price at time t is P{t, T) and which promises to pay one monetary unit at time T. 
In fact, one is left with a residual freedom of choosing A2 [Q. We will return to this 
issue instantly. Let us now restrict the drift Pr to a mean-reverting form: 

p = p{rt) = a- \rt (12) 

Moreover, let us suppose that S(ri,t) = S = const. Let us further assume that 

Tt = r^ + Xt (13) 

Then after transforming from r to x the PDE for the derivative price V reads as: 

(14) 



As pointed out in [|J|, it is not possible to separate the market price of risk A2 
from the difference a = a — A2S. Let us now return to the process for the arbitrage 
return Xt assumed in the last section, Eq.(0). Under the martingale measure for the 
discounted zero bond price, the process for xt obtained from Eq.(|TO]) together with 
(0) and (|13D using standard techniques |25| is given by the SDE: 



dxt ={a- Ar° - \xt) dt + T.dW'^ (15) 

The change to a measure W^ that makes the discounted zero bond price a martingale 
of course amounts to a choice as there is no such instrument in the real world. This 
means that there is no unique martingale measure in terms of real world instruments 
so we are necessarily forced to choose one. One possibility is to require that the 
process Eq. ([T5|) under the martingale measure is mean-reverting to zero: deviations 
from economic equilibrium should disappear to zero. This means that a — \r^ = 0, 
making Eq. ([T5|) identical to the corresponding expression in Eq.(|^). 

Then Eq. (JT^) has a similar but not the same form as Eq. (3) in |Tl| (the last term 
^ instead of +X^ 

ax ax 



should read —Xx^ instead of +A%^ 



1^ + — 5sj + Sg^('-° + -)-(''° + -)'' + Ta^-^- a; = (w) 

and the process for the arbitrage return becomes 

dxt = -Xxtdt + T.dWl (17) 

yielding the consistency of the no-arbitrage approach in the virtual world with the 
arbitrage dynamics proposed in Eq.(|^). The security price dynamics becomes 

dSt = (r° + xt)Stdt + aStdWl (18) 

with respect to the martingale measure. This equation basically couples arbitrage 
returns to security price dynamics under the chosen martingale measure. Thus 
incompleteness is introduced here in terms of stochastic drift as mentioned above. 

Of course, the dynamics of real interest rates are not the same as the dynamics 
proposed here for a virtual world. The relaxational time scale 1/A originating from 
the disappearance of virtual arbitrage returns is much shorter than a time scale of 
mean reversion for real interest rates. As will become clear in the next section, our 
zero bond price in the virtual world will approach a real (constant interest rate) 
bond price in the limit of infinitely fast relaxation dynamics for the arbitrage return 
Xt- One might object at this point, that we have assumed a hedging strategy in the 
virtual world which does not exist in the real world [xt cannot be hedged). Indeed, 
the hedging strategy in the virtual world expressed in terms of the security St and 
a real world cash bond Bf = exp(— r°t) leaves us with an extra amount arising 
from the dynamics of Xt, and is therefore not self- financing in terms of real world 
instruments. We will address this issue in more detail in section 5. 

In order to complete the pricing problem in the virtual world, Eq. ([l^) requires 
the boundary condition, e.g. for a European claim (which is exercised or not exactly 
at time T). It must be chosen as |jTll 



V{t,S,r)U=T = XSix) (19) 



where X is the final payoff (in the real world) depending on S{T) and r = r'^ + x. 
Then the option price in the real world may be calculated as an average over the 
initial arbitrage return as follows: 



V{t, S, r^ 



dxV(t, S, r)p{x) 



(20) 



where p{x) is a probability density function chosen according to the dynamics of Xt 
to be discussed below. 

However, we will not proceed to solve the PDE, but remember that according to 
the Feynman-Kac lemma [|18l O 



V{t,S,r) = E^ 



e-I^dsr,XS{xT)\xt = x,r'^,St 



(21) 



Now it is easy to show (e.g. by using the path integral approach |^) that 



E, 



Q 



e-I^'i^^sxs(xT)\xt = x,r'>,St 



Er 



e-Jld'^^X\xt = x,r^,St]XT 



X p{xT = ^\xt = x,r^,St) 



f22) 



where p{xt = 0|xi = x,r^,St) is the conditional probability density function for 
Xt = given xt = x,r^,St. The last equation allows us to utilize results from 
the literature on option pricing theory with stochastic interest rates. In fact, we 
will solve the pricing problem for the actual payoff function X in an interest rate 
environment where the short rate process in the virtual world r^ starts from r° + x 
at time t and comes back to r" at the time of maturity T, or put otherwise where 
Xt = X and xt = 0. 

The average as given in Eq. (^Oj) will then be performed in a different way. From 
the constraint V{t, S, ro)|t=T = X, it clear that 



p{x) 



p{x) 

m 



(23) 



where p{x) is the probability density function for the initial value x of the arbitrage 
return. Using the fact that in our case p{xt = 0\xt = x,ro, St) = p{xt = 0\xt = x) 

and 

p{x) 



p{xt = x\xt = 0) = p{xt = 0\xt = x) 

p(0) 

one may rewrite the average in Eq.(pOD as follows: 

/oo 
dxV{t, S, r; xt = 0)p{xt = x\xt = 0) 
-oo 

where 



y(t,^,r;xT = 0) = E, 



Q 



/* '^'''X\xt = x,r^,St]XT 







(24) 



(25) 



(26) 



and where p{xt = x\xt = 0) is the conditional probability density for the arbitrage 
return at time t to be equal to x given that its value at T > t, the time of expiry, is 
zero. Its explicit form will be discussed in the next section. In fact, as t = T, one 
obtains p{xt = x\xt = 0) = 6{x) as required. 



4 Valuation of European call and put options in 
the presence of virtual arbitrage opportunities 

Considering the arbitrage return as a part of the stochastic interest rate rj in our 



virtual world, we can now draw upon a classical result of Merton |2g] in order to 
derive formulas for European call and put options. In fact, instead of solving Eq.(|l6|) 
together with Eg. (p!9|) for X = niax(5' — K; 0) (or X = niax(i^ — S"; 0)) for call or 
put options respectively, we consider the security price dynamics Eq.(||) together 
with the following SDE for the price of a zero bond: 

dtP{t, T) = Pit, T) [iipit, T)dt + ap(t, T)dW^) (27) 

Now assuming that o"p(t, T) is a known function of t and T, the price of a European 
call (and put) option c (and p) at time t which expires at time T with strike price 
K is given by [ p^ : 

c = SN{di)-P{t,T)KN{d2) (28) 

p = P{t,T)KN{-d2)-SN{-di) (29) 

where N{x) is the cumulative normal distribution and 

ln(^/i^) - ln(P(t, T)) + a\T - t)/2 



di 



aVT^ 



d2 = di — a\jT — t 
a\T-t) = J ds(a^ + al{s,T)-2psP<yap{s,T)) (30) 

The parameter a is the volatility of the security, and psp is the instantaneous corre- 
lation between the stock and zero bond prices which for simplicity we assume to be 
zero as above. Let us now connect to the stochastic interest rate dynamics rj. If the 
process for P(t, T) is derived from the process for r^ using Ito's lemma, one obtains 
the following dependance of the bond volatility on the parameters of the process for 

n: 

1 dP 

ap(t,T) = ^— --^ (31) 

where S is the short rate volatility. For the specific short rate dynamics chosen in 
Eq.(p!OD and (|12|), the zero bond price can be calculated explicitly as a function of t, 
T, the current short rate level r and the model parameters. Concerning the latter 
ones, the drift of the short rate process has to be risk-adjusted by the market price 
of risk giving a as mentioned above. Let us now calculate the bond price P{t,T). 
Using the fact that 



Pit,T) = E, 



T 



/; dsr. 



' j dsxs 



Xi — X 



(32) 



one obtains from the dynamics of Xt Eq. (p!7| ) (after a tedious calculation given in the 
appendix using the path integral approach pO|) ^he following result: 



P(t,r)=exp(-(r"-^)(T-t)-itanh(^^i^)(^x + ^)) (33) 

10 




Performing the differentiation in Eg . (^Tj) , one obtains: 

ap(t,T) = -|tanh(^^^^) (34) 

Now, tlie restriction a = Xr^ which makes the drift of the process for the arbitrage 
return Xt be equal to —Xxt (under the martingale measure), gives the desired asymp- 
totics of the zero bond price. In fact, as A ^ oo, which can be interpreted as an 
infinitely fast disappearance of virtual arbitrage returns, it reads as 

lim P(t,T) = e-'^"(^-*\ (35) 

A— +00 

the zero bond price for a constant risk- free rate r^. Therefore the restriction on a 
mentioned above, and thus our choice of the martingale measure is reasonable also 
from the viewpoint of correct zero bond price asymptotics. Let us know turn to 
the evaluation of the modified security price volatility a. Evaluating the integral in 
Eg. (pOD , one obtains: 

Likewise, in the limit A -^ 00, the contribution to virtual arbitrage returns disap- 
pears and one recovers the "bare" security price volatility: 

lim a = a (37) 

A^oo 

The asymptotic eguations Eg. (|35| ) and (0) assure that in the case of infinitely fast 
vanishing arbitrage returns the Black-Scholes formulas (for a constant risk-free rate 
r°) are recovered from Eg.(^). When the option approaches maturity, there is the 
following expansion of a^: 

a^ = a^ + ^(T -tf + O ((T - tf) (38) 

Now as the option price in our virtual world is fixed in terms of the parameters of the 
arbitrage return process, we need to turn to the explicit evaluation of the average 
carried out in Eg. (p5[) . For an Ornstein-Uhlenbeck process Xt given by Eg.([T7|) it is 
well known (e.g.[^]), that the transition probability to go from x' at time to a; at 
time t is given by: 



X — x'e^^^ 



p{x, = x\xo = X') = /5 (1 - e-^^Y^' exp I -^I^^—^^L 1 (39) 

What is needed however in our case, is p{xt = x\xt = 0) for T > t which is obtained 
from Eg.([39|) as follows: 

j)(x) 
p{xt = x\xT = 0)=p{xT = 0\xt = x)—-— (40) 

p(0) 
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where p{x) is the probabihty density for x which is obtained as a hmit probabihty 
density form Eq.(B9|) as t — > cxd: 



A X 

^-^exp --. 



P(^) = V^^6xp -— xM (41) 



The final expression for the transition probabihty thus reads: 



P(x. = .|.. = 0) = y'-A- (1 - e--<-'))-"%xp (-^^ \i_,i,.-.,) ) (42) 

It has aU the desired features needed. Using the following representation of Dirac's 
delta function: 

lim rae"""'^' = 5{x) (43) 



one obtains both for the limit of infinitely rapid disappearance of arbitrage returns 

lim p{xt = x\xt = 0) = 6{x) (44) 

A/S2— >oo 

and for the limit t —>■ T oi approaching the option's time of maturity 

lim p{xt = x\xt = 0) = 6{x) (45) 

In both cases, one expects arbitrage returns to disappear. Next, the average over 
virtual arbitrage returns in Eq.(P3[) is carried out explicitly for a European call 
option (for V = c) a.s 

/oo 
dxN{di)p{x{t) = x\x{T) = 0) 
-oo 
/oo 
dxP{t, T)N{d2)p{x{t) = x\x{T) = 0) (46) 

-oo 

and a European put option (for V = p) as 

/oo 
dxP{t,T)N{-d2)p{x{t) = x\x{T) = 0) 
-oo 
^oo 

- S / dxN{-di)p{x{t) = x\x{T) = 0) (47) 



where P{t,T) is given in Eg. (|53D . The integrations with respect to x cannot be 
performed analytically. However, the integrands decrease sufficiently fast to zero as 
X -^ ±cx), so that a numerical integration can be be easily performed. 

It is obvious from intuition that the pricing formulas Eq.s (^) and ( ^Tj) contain 
the fundamental time scale Tarbitrage = 1/-^- In fact, one can introduce the following 
scaled variables: 



^arbitrage 

r 



^^ = A 
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xx = 


X 

A 


crx = 




(T\ = 


a 

Vx 


Sa = 





(48) 

Then A can be eliminated from the pricing formulas. The parameters in Eq. (^0|) can 
expressed in terms of the scaled variables of Eq. 

di 



\n{S/K) - ln(P(u)) + alu/2 



d2 = di - ax\/u 

al = a^ + S^(l-^tanhg)) (49) 



where P{u) is given by: 

P(n) = Pit, T) = exp (- (rl - ^S^) u - tanh (|) {xx + S^ ) (50) 

5 Replicating hedging strategies 

The issue of hedging strategies in the virtual and the real world mentioned above 
will now be addressed. The fact that there is no instrument in the real world to 
hedge intermediate arbitrage returns leads us to conjecture that a hedging strategy 
might not be self-financing. 

To be specific, let us denote a cash bond in our virtual world as follows: 

Bt = exp ( / dsrsj (51) 

It monitors the temporal evolution of the value of an initial cash deposit Bo = 1 
which earns the instantaneous interest rate r^. Let us further introduce the cash 
bond in the real world 

5° = exp(r°t) (52) 

and as a further abbreviation ( which may be termed the "arbitrage bond" ) 

Bf = exp(f dsxs) (53) 

Evidently, one obtains: 

Bt = B^Bt (54) 

Taking Bt for the moment as a real cash bond, a self-financing strategy Vt consists 
of holding (ft in the security St and ipt in the cash bond Bt such that 

Vt = (ftSt + i^tBt => dVt = ifitdSt + iJtdBt (55) 

13 



i.e. the value change dVt is only due to price changes dSt and dBt. For our security 
price model with stochastic interest rates in the virtual world, one can show that 
Eg . (pSD holds |2^. Moreover Vr = X, i.e. the value of portfolio V equals the final 
payoff, i.e. it is replicating. So in terms of our fictitious cash bond Bt there is a self- 
financing, replicating strategy. In the real world, our strategy will remain replicating 
by construction (see Eq.s(^) to (p5|) ). However, it will not be self-financing in terms 
of the real cash bond B^ and the security price St, as can be seen by substituting 



for Bt in Eq.(|| 



dVt = iptdSt + iJtdiB^B^) 

= ^tdSt + AB^dB^ + ^tB^dB^ 

= iftdSt + ^tB^dB^ + ^tB^B^Xtdt 

= ^tdSt + iJtB:dB^ + {Vt-y^tSt)xtdt (56) 

The last step was to replace iptBtB^ = iptBt by Vt — ftSt using Eq. (|55D . The third 
term on the r.h.s of the last line accounts for extra costs or gains due to arbitrage 
opportunities. It is exactly equal to the instantaneous (positive or negative) arbi- 
trage return earned on the delta hedge Vt — ftSt- In fact, one has A = ipt, and 
therefore 

Ut = Vt- iptSt (57) 

where IIj is the delta hedge portfolio discussed in section 2. The replacement r -^ 
r^ + Xt introduced by Ilinski [|10| gives rise to the same additional term in the hedging 



strategy Vt, if one considers the change dllt as follows: 

dUt = (r° + Xt)Utdt = r'^Utdt + xtUtdt (58) 

The second term on the r.h.s of this equation is the source of additional intermediate 
profit and loss (p&l) during the hedging process. Therefore, we conclude that the 
replacement of Ilinski is completely equivalent to the introduction of a fictitious cash 
bond Bt or likewise an interest rate r^ as defined above, which ensures a self-financing 
hedging strategy in the virtual world. 

The additional hedging costs or gains which arise in the real world are covered 
by an additional premium contained in the option price as obtained in Eq. (|25|) (with 
respect to the Black-Scholes price). This premium is positive in most cases as will 
be clarified below when numerical examples are discussed. 

Finally, let us further back up the interpretation of (yt — iptSt)xtdt as representing 
the differential p&l on the hedging strategy within the time interval dt, using the 
following argument (whose formulation is borrowed from [|14|)- At time t, an option 
is sold at Ot in the real world, and using the premium the following portfolio (Vt) is 
bought: 

{Vt) = {ipt)St + {^tB:)B^ = Ot (59) 

where (. . .) corresponds to an average over all paths {xs}s£it,T]- Furthermore, the 
change in wealth of the option seller within the time interval [t, T] in the real world 
is given by: 

AW = Ot + j\^s)dSs + j\Vs - ^sSsVds + j\Vs - VsSs)xsds - X (60) 
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The first term is the option premium earned, the second term gives the cumulative 
gain by the trading the asset, the third one corresponds to the cost/gain of the cash 
bond position (used to finance the position in the asset or set side as excess cash 
respectively) which is proportional to the riskless rate r°, and the fourth term is 
supposed to take into account the p&l due to virtual arbitrage. In fact, the fourth 
term can be added to third term giving an effective cost/gain of the cash bond 
position due to the effective rate r° + x^. Finally the last term is the potential cash 
outflow due to the option's payoff. Now using d{Vt) = {ipt)dSt + {{Vt — iptSt)ij^ + 
Xt))dt one shows that: 



/\W 



Ot + l d{Vs)-X = Ot-{VT)-{Vt)-X = (61) 

as Vt = ^ by construction. As AW vanishes identically, Aiy^(t, 5, r°) (where the 
average is taken as in Eq.(^))) vanishes as well which implies that no intrinsic risk 
remains over the remaining time to maturity of the option, and therefore no risk- 
minimization is necessary. The influence of virtual arbitrage is completely taken 
care of by the option premium. We see also that our hedging strategy in the real 
world is not self-financing at every time step but is self-financing when the time 
integral over remaining life time of the option is taken. 

As discussed e.g. in [^, incomplete markets imply that there is no unique equiva- 
lent martingale measure any more. However, martingale theory may still be used if a 
supplementary constraint is added (see the discussion presented in the introduction) 
which then selects a particular martingale measure. In our case this choice has been 
implicitly made when the arbitrage return becomes part of a fictitious interest rate 
in the virtual world. In fact, both local risk (expected conditional squared cost) and 
replication risk (expected squared deviation of the terminal hedging portfolio to pay- 
off) 1^ are trivially minimized, i.e. zero. A detailed comparison of our approach to 
incomplete markets to the FoUmer-Schweizer approach certainly deserves further 
study. 

6 Some numerical results 

In the following, some results are presented for two market situations, a rather 
incomplete market (FIG.(|I]) and @) and a fairly complete market (FIG.(^ and 
(|)). In the first case, the averaged prices c (and p) , the Black-Scholes prices and 
the payoff functions at maturity are given, for parameters A = 10, T—t = 0.8, S = 2, 
a = 0.2, K = 100, r° = 0.08. The unit of time is 1 year, so A = 10 corresponds 
to the rather long relaxation time Tarbitrage of about 25 trading days, supposing a 
year of 250 trading days. E = 2 is inferred from a daily maximum variation of Xt of 
about 20% in absolute value (at 95% confidence level) according to the discretized 
Langevin equation: 

Ax = xt+i -xt = -XxtAt + XS VAt (62) 

The random variable X is standard normally distributed. Taking At = 1/250, 
X = 1.65 representing the two-sided 95% confidence interval, Xt = (as an initial 
value), one concludes 

E = 9.58Ax (63) 
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For various times to maturity T — t, FIG.(|^) and (D presents differences of c (and p) 
and the Black-Scholes prices for tfie clioice of parameters A = 100, S = 0.4, a = 0.2, 
K = 100, r° = 0.08. A = 100 corresponds to a relaxation time of 2 to 3 trading 
days, whereas S = 0.4 is inferred from a daily variation of Xf of about 4% in absolute 
value. 

Let us now comment on the results. Focusing first on the qualitative behavior, 
over a reasonable range of the moneyness parameter m = S/K, the price of a 
European call or put option (c or p respectively) under the influence of virtual 
arbitrage is higher than the Black-Scholes value (see FIG.(||) and (^). The difference 
is more pronounced at the point of maximum curvature which is around m ~ 1 or 
below , whereas it decreases whenever m < 1 or m > 1 (see FIG.(^ and (^). For 
m ^ 1, the call option price is less than the Black-Scholes value. As the time 
to expiry increases the positive difference (except for m ^ 1) increases, and the 
maximum difference is shifted to lower values of m. 

Leaving aside for the moment the negative difference appearing for call option at 
m ^ 1, it appears reasonable that the existence of virtual arbitrage returns causes 
the option price to be above the Black-Scholes value, as deviations from equilibrium 
in general lead to an increase in hedging costs, i.e. the costs for readjusting a 
replication portfolio which is supposed to provide for the final payoff of the option. 
This effect needs to be accounted for in the option premium. The fact that the 
absolute difference to the Black-Scholes result is the largest at the point of maximum 
curvature of the pricing function is understandable from the F ( "gamma" ) risk point 
of view. F denotes the second derivative of the option price with respect the asset 
price S and gives a measure for the non-linear dependance of the option on the 
underlying asset. This non-linear risk inherent to options can be only be hedged by 
buying or selling other options. Any deviations from financial equilibrium due to 
arbitrage opportunities will affect both the option at hand and the options chosen for 
hedging. Moreover, the additional term arising in Eq.(^) leading to intermediate 
P&L during the hedging process is proportional to the delta hedge Vt — iptSt which 
is most relevant at the point of maximum F where the delta hedge is insufficient. 
Therefore, these numerical results are completely consistent with our mathematical 
discussion of the hedging strategy. 

The influence of intermediate arbitrage returns grows as the time to expiry of the 
option increases on the scale of Tarbitrage (see FIG.(|^) and (|D) (all other parameters 
being constant). Several deviations from equilibrium during the life time of an option 
seem to accumulate leading to a higher additional risk premium on the option price. 

Returning to the issue of the negative difference for call options that are far 
in the money m ^ 1 in FIG.(|1|) and (^), a possible explanation is an "overheated" 
market, where deviations from equilibrium tend to relax from the current asset price 
to a lower equilibrium price. This information is accounted for by pricing the option 
at a discount with respect to the Black-Scholes value at the current asset price: the 
market is expected to decrease to a lower price level. 

Considering the quantitative differences between c (and p) and the Black-Scholes 
prices for calls and puts, they are obviously more pronounced in an imcomplete 
market (FIG.(0) and (^), than in a rather complete market where arbitrage returns 
are small and relax fast ((^ and (^)). The numerical analysis given here may be 
refined in various ways (according to the parameter dependances of the options 
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prices) which is the subject of future work. 

As opposed to our resuhs, e.g. the first order correction to the Black-Scholes 



prices for calls given in [10| increases monotonously with moneyness m. Obviously, 



as our result is based on the same model as in []T0| (see section 5) and is exact (apart 
from the remaining integration over the initial arbitrage), some error is made in the 
perturbative treatment. Let us point out here again that the Ornstein-Uhlenbeck 
dynamics of arbitrage returns used here and the fact that the Xt gives the extra 
return on the delta hedge Vt — ftSt, it is quite reasonable that the difference to the 
Black-Scholes price should show a maximum at the price level where the delta hedge 
fails. 



7 Conclusion 

Using the (Ornstein-Uhlenbeck type) relaxational dynamics for "virtual" arbitrage 
returns introduced in [|10[, we have derived closed formulas for simple ("plain vanilla") 
European calls and puts in the presence of arbitrage opportunities appearing and 
disappearing on an intermediate time scale Tarbitrage = 1/A. This result which has 
not been derived previously is obtained using martingale option pricing theory for 
incomplete markets (in the sense of 0]), by making the arbitrage return process 
part of an interest rate process in a virtual world. The influence on option prices 
in the real world (in the presence of rapidly appearing and disappearing arbitrage 
opportunities) is taken into account by summing over the initial arbitrage return, 
and imposing the constraint that arbitrage is absent at the time of maturity of the 
option. 



Comparing our work to |10| , |TT|, flrst, we consider the analysis given above as 
conceptually more clear as to where arbitrage-free pricing fails and where it does not. 
Therefore, in the present work a different route has been proposed by introducing 
a second source of randomness in the derivative pricing problem (apart from the 
security S) right from the beginning. As a consequence, a two variable version of 
Ito's lemma must be used, giving a PDE equation for the derivative price in a virtual 
world which is flnally summed over xt to yield the real world price. Second, instead 
of making the constraint that arbitrage return should vanish at maturity a part of 
the payoff function in the virtual world as in |]10| , we enforce it when the average over 



virtual arbitrage return is taken. This procedure allows us to proflt from Merton's 



classical result on option pricing in a stochastic interest rate environment |^ and to 
arrive at closed-form (up to a numerical integration over the initial arbitrage return 
which is easy to perform) pricing formulas for simple European call and put options. 
Furthermore it has been shown that any hedging strategy will not be self- 
flnancing in the real world where the arbitrage return is not directly observable. 
However, on the average any intermediate costs arising during the hedging process 
are covered by an additional premium contained in the option price. In this sense, 
a hedging strategy can be found that is self-flnancing in a time average sense, i.e. 
when summed over the remaing life-time of the option. The derivation of pricing 
formulas rests crucially on the selection of a speciflc measure from a set of equiva- 
lent martingale measures that contains more than one element, due to intermediate 
market incompleteness which arises because of virtual arbitrage opportunities. 
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The present work may be extended in various directions. The relaxational dy- 
namics of the arbitrage return may be considered to be more comphcated as proposed 
here where it follows a simple Ornstein-Uhlenbeck process. However, additional 
model parameters introduce more sources of model error from the practitioner's 
point of view as each parameter has to calibrated to the market. Furthermore, the 
constraint xt = 0, i.e. that arbitage returns should disappear at the time of ma- 
turity of the option, may be relaxed to allow for a hedging mismatch at maturity. 
This amounts to give up the constraint that the hedging strategy is replicating. The 
extension of this work to the case of correlations between the asset price St and the 
arbitrage return Xt is under way. Certainly, the comparison of the present model to 
stochastic volatility models deserves further study. 
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Appendix 

Following [^] we propose to evaluate the expectation value 



I = Er 



e -It "ixf = X 



(64) 



The expectation value can be stated in terms of a quotient of path integrals as 
follows: 



V ^ 1 L = _ (55) 



Now the numerator and the denominator can be mapped to the propagator of the 
harmonic oscillator in the presence of an external field, and can thus be evaluated 



28[]. The expression for the numerator reads as 



^ Vcosh(A(r-t)) 



where 



2S2sinh(A(T-t)) 
+ 2(e^(^-*)-l)(Cx + C2)) 

+ ^^'1 (66) 



C = ^ (e-(-*) - 1) (67) 
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Likewise one obtains an expression for the denominator: 



^ = V 27rS^ sinh(A(T - t) ^"^ (2^"' " 2S^ sinh(A(T - t)f ^"^^^^^^^ " '^^) 
Calculating X/Y gives the result 






^ = exp(;Ti^(r-t)--tanh^— ^ x+- (69) 



which leads to Eg. (|33D . 
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FIG. |l]: The call option price as a function of moneyness m = S/K (dashed curved 
line: with virtual arbitrage; solid line: Black-Scholes formula). The dashed straight 
line is the payoff function at maturity. Parameters: A = 10, T — t = 0.8, S = 2, 
a = 0.2, K = 100, r° = 0.08. 

FIG. ^: The put option price as a function of moneyness m = S/K (dashed curved 
line: with virtual arbitrage; solid line: Black-Scholes formula). The dashed straight 
line is the payoff function at maturity. Parameters: see FIG. [l|. 

FIG. 1^: Difference of the call price to Black-Scholes value (in absolute value) for var- 
ious values T-t. Other parameters: A = 100, S = 0.4, a = 0.2, K = 100, r° = 0.08. 

FIG. ^: Difference of the put price to Black-Scholes value (in absolute value) for 
various values T — t. Other parameters: see FIG. |]. 
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